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ABSTRACT 

We have studied a sample of Large Magellanic Cloud red giant binaries that lie 
on sequence E in the period-luminosity plane. We show that their combined light and 
velocity curves unambiguously demonstrate that they are binaries showing ellipsoidal 
variability. By comparing the phased light and velocity curves of both sequence D and 
E variables, we show that the sequence D variation - the Long Secondary Period - is 
not caused by ellipsoidal variability. We also demonstrate several further differences 
between stars on sequences D and E. These include differences in velocity amplitude, 
in the distribution of eccentricity, and in the correlations of velocity amplitude with 
luminosity and period. We also show that the sequence E stars, unlike stars on sequence 
D, do not show any evidence of a mid-infrared excess that would indicate circumstellar 
dust. 
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1 INTRODUCTION 

Long Period Variables (LPVs) are known to fall on different 
sequences in the period- luminosity plan e (Wood et al. 1999 ; 



Soszyriski et al . 2004a; llta et al 



ISoszvnski et alj [20071 ; iFraser et al 



y plan e ( Wood et al. i^TO ; 
I [20041; IFraser et alj liooj 



2008). One of these se- 



quences, known as sequence E, is thought to consist of red 
giants in close binary systems showing ellipsoidal variability, 
although this has not been unambiguously demonstrated. A 
small number of the binaries appear to be eclipsing and oth- 
ers appear to have unexpectedly eccentric orbits, based on 
their light curve shape (Soszyriski et al. 2004b). 

A star in a close binary will have its Roche Lobe dis- 
torted by the tidal influence of its orbiting companion. When 
the star begins to fill its Roche Lobe it takes on an elon- 
gated, or ellipsoidal shape, becoming an ellipsoidal variable. 
The velocity variations of such stars are dominated by the 
orbital motion, but light variability is caused mainly by the 
change in apparent surface area as the star orbits around its 
companion. Because of this, the light curve of an ellipsoidal 
variable shows two maxima and minima per orbit - two cy- 
cles for every one cycle of the velocity curve. This is an easy 
way to unambiguously identify ellipsoidal variables and we 
use this test here. 
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Four of the other sequences of LPVs - A, B, C' and C 

- are known to harbour radially pulsating variables (Wood 
et al. 1999). A fifth sequence, sequence D, contains stars 
which show Long Secondary Periods. The origin of Long 
Secondary Periods (LSPs) remains something of a mys- 
tery, though sever al attempts have been made to discover it 
(Wood et al 1999 ; [Hinkle et al.ll2002l; lOlivier fc Woodll2003l ; 
IWood et all |2004| ; iNicholls et al.ll2009h . Due to their over- 
lap in the period-lum inosity diagram , some authors (e.g. 
Soszyhski et al. 2004b, ISoszvriskill2007h have suggested that 
stars on sequences D and E may be fundamentally the same 

- binaries showing ellipsoidal variability. 

Here we study a sample of sequence E binaries. In par- 
ticular we present new radial velocity data derived from 
VLT spectra, which we use alongside MACHO light curves. 
We show that the variations of stars on sequences D and E 
are caused by different mechanisms. Prelim inary results for 
three stars are given in lAdams et al.l |2006). 



2 OBSERVATIONS AND DATA REDUCTION 

The observations and data reduction for this sample are 
the same as for our rec ently published sequence D sam- 
ple (|Nicholls et al.1 |2( )()!) :■ . The spectra were taken using 
the FLAMES/ GIRAFFE spectrograph (Pasquini et al. 
2002) on the European Southern Observatory's Very Large 
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Table 1. Radial Velocities of Sequence E Stars. Stars are identi- 
fied by their MACHO numbers. 



our s equence E red giants (see section 2.5 of iNicholls et al.l 
2009). We calculated the mass function, 



HJD 


77.7429.189 


77.7548.68 


77.7672.98 


77.7673.79 


2954.8508 


0.00 


247.1836 


270.6885 


265.7394 


3005.8604 


0.00 


242.5239 


260.0903 


258.5819 


3067.6028 


0.00 


237.8895 


225.7872 


239.5893 


3091.5571 


0.00 


238.4238 


251.1337 


265.9271 


3280.8611 


272.6101 


255.7782 


284.0382 


244.0299 
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Figure 3. A histogram of the velocity amplitude for our sample. 



Telescope (VLT), on 21 nights from 2003 November to 
2006 March. Radial velocities were calculated via cross- 
co rrelation with t he ir af task fxcor. The reader is referred 
to lNicholls et all l|2009h for details. TableHJshows the radial 
velocities calculated for part of our sample for a few dates. 
The full table, with radial velocities of our whole sample for 
all dates, is available online. 



3 RESULTS 

Plots of phased light and velocity variations for all the stars 
in our sequence E sample are shown in Figs.[T]and[2] A clear 
doubling of the phased light curve with respect to the phased 
velocity curve can be seen for all stars in our sample. This 
behaviour unambiguously demonstrates that these stars are 
ellipsoidal variables. Many of the stars show lightcurves with 
equal maxima and unequal minima, which suggests that one 
end of the ellipsoid is hotter and brighter than the other. 
From Figs. [T] and [2] we see that the deepest minimum oc- 
curs at mean velocity during decreasing radial velocity, i.e. 
when the red giant is behind its companion. Therefore the 
inner end of the red giant ellipsoid - the end closest to the 
companion - is the cooler end, and the side furthest from 
the companion is hotter. This is most likely due to gravity 
darkening towards the companion. 

A histogram of the velocity amplitudes of our sample is 
shown in Fig. [3] The majority of values lie between 15 and 
55 kms -1 , and the mean velocity amplitude is 43.3 kms -1 . 

We have made a binary fit to the velocity curves of 



f(m) = 



K P 

2tvG 



{m + M) 2 ' 



(1) 



using the observed values for period (P) and velocity semi- 
amplitude (K) (Table[2)| . Using the calculated mass function 
of our binary fit and assuming a total system mass of 2 Mq , 
we calculated an estimated companion mass for each of the 
stars in our sample. As expected, for most systems the com- 
panion is less massive than the red giant, but for two stars 
this resulted in companions with significantly higher mass 
than the red giant. This is unlikely since the more massive 
companion should evolve to the red giant stage first, and 
even if the current red giant was originally the less massive 
star, its companion would most likely now be a white dwarf 
of lower mass than the red giant. In these two cases we chose 
the total mass to be such that the red giant and its compan- 
ion had equal mass. This results in higher masses for these 
stars (total system masses of ~ 4.5Mq.) The mass estimates 
are shown in Table [5] alongside the other parameters of the 
binary fit to the velocity. Here 7 is the system velocity, e is 
the eccentricity, to is the angle of periastron, T is the date 
of periastron, asini is the semimajor axis of the red giant's 
orbit, f(m) is the mass function, M is the mass of the red 
giant, and m the mass of its companion. We assume sin i = 1 
when calculating the masses. The errors for each element are 
shown on the line below each star. 

The eccentricity of a binary system is exp ected to de - 
crease over time due to circularising tidal forces (|Zahnll977M . 
provided no mechanism is in place that will stabilise or in- 
crease the orbital eccentricity. Thus binaries with one or 
more red giant components are expected to describe rela- 
tively circular orbits. We confirmed this theory by calculat- 
ing the c ircular i sation time for our stars, using the formula 
given m ISokerl <|200Ch . The median circularisation time of 
our sample is only ~ 3500 y. For comparison, the time a 
1 Mq, Z = 0.008 red giant takes to double its radius when 
R ~ 30 Rpi is ~ 2 x 10 7 y (using the evolutionary tracks 
of lGirardi etaT] 2000). The value of 30 R© was chosen as a 
typical Roche Lobe radius for these stars (see the velocity 
amplitude calculations given later in this section). A dou- 
bling of the radius is an estimate of the evolution time over 
which the tides have had time to act. Ten of the eleven stars 
in our sample have e < 0.1, and the mean eccentricity of 
the sample is 0.07. Therefore the majority of our sample of 
red giant binaries have fairly circular orbits, as expected. 
However some sequence E stars do show confoundingly high 
eccentricities (see Soszyhski et al. 2004b), and this is a prob- 
lem that must be solved in the future. 

The relation between velocity amplitude and period for 
our sa mple and for the sequence D sample of INicholls et al.l 
(2009) is plotted in Fig. [4] In a binary system, for given 
masses, a longer period means a wider orbit, and thus a 
smaller amplitude in velocity variation. Fig.|4]shows velocity 
amplitude decreases with increasing period for the sequence 
E stars, as expected. However the sequence D stars show the 
same velocity amplitudes for all periods and do not follow 
the expected binary relation. A vivid demonstration of the 
disparity between these two samples is found in the region 
where their periods overlap. At periods between 200 and 
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Figure 1. Phased MACHO blue (Mb) and red (Mr) light curves and phased radial velocity curves for our sequence E stars. The blue 
line shows the binary fit to the velocity curve. 
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Figure 2. Same as Fig. [T] 
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Table 2. Orbital Elements for Sequence E Stars 



Star 7 K e u T Pa sin i f(m) m 

(kms" 1 ) (kms- 1 ) (deg) (HJD) (days) (R ) (M ) (M~) (M ) 



77.7429 


189 


281.40 


25.99 




0.05 


275 


38 


3332.5 


110 


34 


56.58 







2003 


1 


07 





93 






± 0.57 


± 0.80 


± 


0.03 


± 28 


52 


± 8.8 






± 1.75 


± 





0186 










77.7548 


68 


247.20 


8.54 




0.07 


170 


67 


3502.6 


442 


83 


74.55 







0285 


1 


52 





48 






± 0.19 


± 0.24 


± 


0.03 


± 25 


63 


± 32.9 






± 2.22 


± 





0024 










77.7672 


98 


250.71 


32.77 




0.06 


38 


99 


2989.4 


156 


13 


100.92 







5679 


2 


27 


2 


27 






± 0.75 


± 1.10 


± 


0.03 


± 29 


11 


± 12.3 






± 3.42 


± 





0575 










77.7673 


79 


252.37 


24.01 




0.02 


236 


19 


3064.4 


131 


35 


62.30 







1887 


1 


09 





91 






± 0.41 


± 0.59 


± 


0.02 


± 61 


77 


± 22.6 






± 1.55 


± 





0140 










77.7789 


152 


255.04 


14.03 




0.05 


198 


55 


3024.2 


157 


73 


43.67 







0451 


1 


41 





56 






± 0.74 


± 0.97 


± 


0.08 


± 90 


81 


± 40.5 






± 3.05 


± 





0094 










77.7790 


72 


275.00 


15.35 




0.19 


262 


20 


3326.7 


328 


79 


97.85 







1167 


1 


22 





78 






± 0.83 


± 0.92 


± 


0.06 


± 13 


11 


± 10.2 






± 6.02 


± 





0214 










77.7791 


115 


350.88 


21.47 




0.03 


275 


04 


3032.0 


139 


18 


59.02 







1429 


1 


17 





83 






± 0.30 


± 0.35 


± 


0.02 


± 40 


51 


± 16.3 






± 0.97 


± 





0070 










77.7910 


41 


238.06 


20.41 




0.08 


10 


03 


3075.8 


191 


46 


76.94 







1673 


1 


13 





87 






± 0.55 


± 0.66 


± 


0.04 


± 21 


82 


± 11.8 






± 2.53 


± 





0163 










77.7910 


77 


255.34 


44.35 




0.08 


279 


65 


2987.6 


65 


72 


57.38 







5890 


2 


36 


2 


36 






± 1.67 


± 2.83 


± 


0.04 


± 38 


55 


± 7.1 






± 3.67 


± 





1131 










77.7912 


111 


231.16 


14.76 




0.06 


223 


73 


3075.1 


128 


89 


37.50 







0428 


1 


11 





56 






± 0.45 


± 0.62 


± 


0.04 


± 41 


73 


± 14.7 






± 1.57 


± 





0054 










77.7914 


74 


305.31 


16.33 




0.02 


272 


80 


3271.8 


282 


24 


91.04 







1276 


1 


20 





80 






± 0.11 


± 0.14 


± 


0.01 


± 24 


26 


± 18.9 






± 0.79 


± 





0033 











a Here we assume sin i = 1 . 
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Figure 4. Full velocity amplitude plotted against orbital period 
for ou r sequence E sample (blue points) and the sequence D sam- 
ple of lNicholls et alj |2009) (red crosses). 



450 days, the difference in velocity amplitude between the 
sequence D and sequence E samples is as much as 30 kms -1 . 

Ellipsoidal variability is only visible once a star has sub- 
stantially filled its Roche Lobe. Stars in wider orbits (with 
lower velocity amplitudes) will fill their Roche Lobes when 
they are further up the giant branch, and thus will be more 
luminous than stars in closer orbits. For a given luminosity 
(and radius) there is a range of velocity amplitudes the star 
may have, the maximum of which is dictated by the size of 



the closest possible orbit the companion can occupy. There- 
fore, for our ellipsoidal variables we expect that the max- 
imum velocity amplitude should decrease with increasing 
luminosity, but that stars may occupy velocity amplitudes 
below the maximum for a given luminosity. This is shown 
in Fig.[5]which gives a plot of K magnitude against velocity 
amplitude for our sequence E stars. In order to define the 
upper limit of velocity amplitude for a given luminosity, we 
calculated the minimum orbital separation for a theoretical 
sample of binaries with equal mass components (mass ratios 
of q = 1) in which the red gian t is filling its R oche Lobe. 
Using the approximation given in Egglcton (1983), we calcu- 
lated the Roche Lobe radius tl to be 0.37, in units of orbital 
separation. Therefore the minimum orbital separation for a 
Roche-Lobe filling binary with q = 1 is a — = 2.7tl - We 
substitute for rz, the radius expected for a given luminosity 
in these stars, calculated from a fit made to the radius-if 
magnitude data for our sample. Finally we calculated the 
maximum velocity amplitude assuming a circular orbit and 
components of mass 1 M . This velocity amplitude upper 
limit is shown by the solid blue line in Fig. [5] 

Most stars lie where expected in Fig. [5] but two stars 
lie above the upper limit line. These are the two systems 
with higher mass components, mentioned earlier. As both 
these systems have components of around 2.3 M , we also 
calculated the maximum velocity amplitude for stars of this 
mass, and this is shown by the green dashed line in Fig. [5] 
One star lies above this limit. It appears somewhat atypical, 
as it has a significantly higher effective temperature than 
the rest of the sample (~ 5200if compared to the median 
4200^) and hence a smaller radius. It therefore does not 
overflow its Roche Lobe as Fig. [5] suggests. 
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Table 3. Sequence E Star Properties 





Star 


R 


L 


Teff 


a 






(R©) 


(Lq) 


(K) 


(Re) 


77 
1 I 


1 A on ICQ 

/ 4zy. ley 


a c no 


oyo.it) 


A(\*7~\ A Q 

4Uf 1.4o 


IZl.oo 


77 


7548.68 


92.03 


1917.72 


3980.39 


307.67 


77 


7672.98 


56.13 


869.34 


4181.83 


201.84 


77 


7673.79 


42.24 


591.59 


4378.57 


136.84 


77 


7789.152 


52.85 


841.10 


4274.26 


154.60 


77 


7790.72 


92.71 


1736.69 


3868.59 


252.27 


77 


7791.115 


49.66 


805.48 


4362.15 


142.23 


77 


7910.41 


72.37 


1790.14 


4412.02 


175.91 


77 


7910.77 


38.58 


983.43 


5202.34 


114.76 


77 


7912.111 


47.41 


604.76 


4155.77 


135.12 


77 


7914.74 


70.70 


1079.15 


3933.14 


227.86 
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Figure 5. K magnitude (from 2MASS) plotted against velocity 
amplitude for our sequence E stars. The point size shows relative 
light amplitude in MACHO red. The solid blue curve gives the 
upper velocity limit for equal- mass components of 1 Mq, and the 
dashed green curve gives the upper limit for equal-mass compo- 
nents of 2.3 Mq. See text for details. 

Table [3] gives the radius, luminosity, effective temper- 
ature and orbital separation data for our sample. For the 
majority of stars, these properties were calculated using 
the 2MASS J and K magnitudes. However for the star 
77.7910.77, which lies outside the green dashed limit in 
Fig. [S] the properties were calculated from OGLE V — I val- 
ues, which are more reliable for warmer stars. Due to this 
star's warmer temperature, we suspect it could be a helium 
core burning clump star or a blue loop star. 

3.1 The mid-infrared colour of sequence E stars 

It has recently been shown (|Wood fc Nicholls! 12009) that 
variable red giants belonging to sequence D have a mid- 
infrared excess when compared to similar red giants with- 
out the Long Secondary Period that characterises sequence 
D stars. This mid-infrared excess is assumed to arise from 
circumstellar dust associated the presence of the LSP. Since 
the sequence E stars are close binary systems, some of which 



have quite eccentric orbits (Soszyhski et al. 2004b), there is 
a possibility that these systems could have substantial cir- 
cumstellar disks of dust and gas. 

In order to investigate this possibility, we searched for 
a mid-infrared excess in sequence E stars in the LMC using 
mid-infrared data from the Spitzer Space Telescope SAGE 
survey (Blum et al. 2006; Meixner et al. 2006). The se- 
quence E stars we re obtained from the catalogue given by 
iFraser et al] (|2008l h None of these objects were detected at 
24 fim in the SAGE survey but 262 were detected at 8 /im. 
The light curves of these 262 stars were examined and it 
was found that only 184 could be considered as definite se- 
quence E stars, i.e. ellipsoidal or eclipsing binary systems. 
The remainder were variables whose light curve character- 
istics showed that they belonged to sequence D or the pul- 
satio n sequences 1 to 4 (using the notation in IFraser et al.l 
l2008h or they were stars of indeterminate light curve type. 

In order to see if the sequence E stars show a mid- 
infrared excess, we obtained a comparator sample of stars 
that were similar to the sequence E stars, but without binary 
companions. Such stars are the normal, non-varying field 
red giants in the LMC. Since the sequence E stars that were 
detected at 8 /im have 12.5 < K < 14 and 0.6 < J - K < 
1.4, we selected comparison sources from the SAGE survey 
with these characteristics along with the requirement of an 
8 fim detection. An examination of this sample of ~40000 
objects showed that sources fainter than K = 13.5 were 
near the 8 /im faint detection limit and that this significantly 
biased the population of detected stars in favour of objects 
with a higher K-[8] colour. We therefore omitted such stars 
and we compared sequence E stars and non- varying field red 
giants only in the interval 12.5 < K < 13.5 . Fina lly, all 
variable stars in the catalogue of IFraser et al.l (|2008f i were 
removed from the list of comparison sources. This left 144 
sequence E stars and ^33000 non-varying field red giants. 
The K-[8] colour distributions of these two samples were 
then examined to search for an 8 fim mid- infrared excess. A 
two sample K-S test gives a probability of up to 0.89 that 
the sequence E stars and non-varying field red giants come 
from the same underlying distribution in K- [8] colour. In 
other words, there is no evidence that sequence E variability 
generally leads to a significant amount of circumstellar dust 
and a mid-infrared excess. 



4 DISCUSSION AND CONCLUSIONS 

With new velocity curves, we now have strong evidence that 
the sequence E variation is indeed caused by ellipsoidal vari- 
ability. This is shown by the doubling of the light curve with 
respect to the velocity curve in these stars, as in Figs. [T] 
and[2] 

However, as we showed in iNicholls et ail (|2009h , the 
phased light curves of sequence D variables - the stars with 
Long Secondary Periods - do not show this doubling phe- 
nomenon (see fig. 3 of that paper) . Although it has been sug- 
gested that the sequence D and E variations may have a com- 
mon origin due to their proximity in the period-luminosity 
plane (Soszyhski et al. 2004b), and the possible existence 
of ell ipsoidal shapes in their residual light curves ( Soszyhski 
l2007h . no sequence D star has so far been found whose light 
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executes two cycles during one velocity cycle. This seems to 
rule out ellipsoidal variation in sequence D stars. 

A further difference between the stars on sequences D 
and E is demonstrated by their respective velocity ampli- 
tudes. The sequence E stars of the current sample show 
full velocity amplitudes of at least 15 kms - some much 
larger - as can be se en in Fig. [3] However as we showed in 
iNicholls et all (l2009). stars with LSPs have much smaller ve- 
locity amplitudes, typically around 3.5 kms -1 . This marked 
difference in the size of the velocity variations between se- 
quence D and E - which can be seen in Fig. [4] - supports 
our assertion that they are caused by different mechanisms. 

The sequence E stars, known binaries, show an expected 
spread in companion mass estimates. However this distribu- 
tion seems very different for the sequence D stars. If treated 
as binaries, sequence D stars have companions that are small 
and absurdly similar (~ 0.09 Mq). Additionally, the angles 
of periastron of sequence D stars are heavily biased towards 
large values, whereas we would expect a uniform distribution 
for binaries. Unfortunately we cannot compare the distribu- 
tion of angle of periastron for the current sequence E sample: 
as most of these stars have almost circular orbits, angle of 
periastron is poorly defined (see the large errors in Table [2]). 

The distribution of eccentricity is very different for se- 
quence E and D stars, when the latter are treated as bina- 
ries. Most of the red giant binaries on sequence E have low- 
eccentric ity orbits, as e x pecte d for their evolutionary state. 
However INicholls et alj (|2009r ) showed that the sequence D 
stars have a significantly higher eccentricity, if we assume 
they are binaries (the mean eccentricity of that sample is 
0.3). A two-sample K-S test gives a probability of less than 
2.1 x 10 -6 that the sequence D and E eccentricity distribu- 
tions come from the same underlying distribution. 

Further evidence that sequence D stars are not el- 
lipsoidal variables is given by the relation between lumi- 
nosity and velocity amplitude in those stars. In Fig. [5l 
we showed that the sequence E stars generally populate 
an area delineated by a maximum velocity amplitude- 
luminosit y relation, as expec ted for ellipsoidal variables. 
However. INicholls et al.l (|2009l ) showed that the sequence D 
stars do not show any correlation between these properties. 

Finally, we have searched for a mid-infrared excess in 
sequence E stars since such an excess was found among the 
sequence D stars. No mid-infrared excess was found, demon- 
strating another difference from the sequence D stars. 

In summary, we have demonstrated several significant 
differences between stars on sequences D and E. Most par- 
ticularly, we have shown that Long Secondary Periods - the 
sequence D variation - are not caused by ellipsoidal variabil- 
ity, unlike the variation of the sequence E stars. The cause 
of Long Secondary Periods remains, at this time, a mystery. 
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